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Abstract 

A preliminary flux estimate of various cosmic-ray constituents based on the atmospheric 
Cerenkov light flux of extensive air showers using fractal and wavelet analysis approach 
is proposed. Using a Monte-Carlo simulated database of Cerenkov images recorded by the 
TACTIC telescope, we show that one of the wavelet parameters (wavelet dimension Bq) 
provides > 90% segregation of the simulated events in terms of the primary mass. We use 
these results to get a preliminary estimate of primary flux for various cosmic-ray primaries 
above 5 TeV energy. The simulation based flux estimates of the primary mass as recorded 
by the TACTIC telescope are in good agreement with the experimentally determined values. 
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1 Introduction 



One of the main problems in the field of primary cosmic-ray investigations is the 
form of the energy spectrum and the precise determination of the mass compo- 
sition. Precise information about the energy spectrum and mass composition of 
primary cosmic-rays at VHE (> 10 11 eV) energies is an important input for un- 
derstanding the origin of cosmic -rays and their propagation through the Galaxy. 
Cosmic -rays exhibit a nearly smooth power law energy spectrum over more than 
10 decades of energy (from ~ 10 9 to 10 20 eV) with significant kinks at ~ 10 15 eV 
(knee) and ~ 10 18 eV (ankle)[l,2]. The knee feature at ~ 10 15 eV, which involves 
a change in the differential spectral exponent from ~ 2.7 to ~ 3.1 has defied a 
clear explanation so far. Among the possibilities, one of the reasons suggested for 
the knee feature involves a change in the composition of the primary cosmic-rays 
from a mainly proton dominated(~ 90% protons) to a mixed (protons + heavy 
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nuclei) composition or even heavy nuclei domination[2,3]. Thus, it becomes very 
important to measure the cosmic-ray composition in the VHE range, especially at 
energies around 10 14 -10 15 eV. Although primary cosmic-ray composition is known 
more or less precisely at energies up to ~ 10 TeV from direct satellite/balloon based 
observations[4,5], the low flux of primary cosmic -ray particles at energies > 10 13 
eV (~ 10 _2 particles cm~ 2 s _1 )does not allow direct measurement of primary com- 
position with detectors carried on balloons and satellites. Attempts have, therefore, 
been made to infer the cosmic-ray composition indirectly from data collected with 
large air shower arrays and atmospheric Cerenkov telescopes which detect primary 
cosmic -rays indirectly by recording the secondary particles and photons produced 
in the atmospheric cascades initiated by these primaries[6,7,8,9]. These attempts 
have achieved limited success so far because of the difficulty in relating the ob- 
served cascade parameters to the primary mass. In the present communication, we 
make an attempt to show that one of the fractal/wavelet parameters derived from 
the Cerenkov images recorded by an imaging Cerenkov telescope (like TACTIC) 
exhibits an unprecedented efficiency for separating primary cosmic-rays in terms 
of their mass. This makes it possible to estimate the flux of cosmic-ray primaries at 
energies > 10 13 eV at large zenith angles with ground-based atmospheric Cerenkov 
telescope systems. 



2 Imaging Atmospheric Cerenkov Technique 

Primary cosmic -rays (protons, hadrons, photons etc.) with energy > 10 11 eV, initi- 
ate extensive air showers in the atmosphere, generating in the process a large num- 
ber of secondary particles (e ± J fi ± 7 n ± , hadrons, photons) some of which survive 
absorption in the atmosphere and can be detected at ground level. A large number of 
e ± , traveling with relativistic speeds, generate Cerenkov photons in the atmosphere 
which arrive in the form of a few ns duration optical pulse (A~ 300 — 600nm) 
spread over a radius of ~ 100m at ground level[6,7,10]. An imaging Cerenkov 
telescope, comprising a large light collector with a matrix of photomultiplier tube 
(PMT) detectors in its focal plane, generates a two-dimensional picture of the re- 
sultant Cerenkov photon distribution (called the Cerenkov image ), the geometrical 
features of which depend upon the nature (mass) of the cosmic-ray primary which 
initiated the cascade. The technique has been successfully employed by various 
imaging Cerenkov telescopes including TACTIC to detect VHE gamma-ray signals 
from a number of galactic and extragalactic sources[7,8]. The TACTIC (TeV Atmo- 
spheric Cerenkov Telescope with Imaging Camera) telescope at Mt.Abu (24.6° N; 
72.0° E; 1300m asl) comprises a 3.6m aperture, tessellated light collector of overall 
parabolic shape mounted on an altitude- azimuth mount and equipped with a 349- 
pixel imaging camera at its focal plane [10,11]. The imaging camera provides a field 
of view of 3° radius with a uniform pixel resolution of 0.31°. Using a topological 
trigger based on a coincidence^ ~ 20ns) between nearest-neighbour, non-collinear 
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triplets (3NCT) in the central 225 pixels (15x15), the imaging telescope records 
Cerenkov images of extensive air showers initiated by primary gamma-rays of > 
1.5 TeV and primary hadrons of energy > 3 TeV from the vertical direction [10]. 
The images are composed of Cerenkov photon intensity observed by each pixel in 
the form of charge to digital converter counts(CDC). The effective collection area 
(area of Cerenkov light pool at ground level) of the system has been estimated to be 
~ 3 x 10 8 cm 2 [17]. Both the threshold detection energy and the effective detection 
area increase with the zenith angle of observation[10,12]. 



3 Simulations 

The results presented in this paper are based on the analysis of simulated Cerenkov 
images produced by several primary species like gamma-rays, protons, Ne and Fe 
nuclei. The simulations have been performed for the 349— pixel imaging camera 
of the TACTIC telescope using the CORSIKA (version 6.51) air shower simula- 
tion code[14,15,16], supplemented by appropriate custom built software to take 
into account atmospheric absorption, reflection at the TACTIC light collector sur- 
face, PMT characteristics and the night sky background induced shot-noise fluctu- 
ations[7,13]. For each primary species 6,00,000 showers in different energy ranges 
( 2TeV, 4TeV, 5TeV, 8TeV, lOTeV, 15TeV, 20TeV, 30TeV) with energy spectral in- 
dex of -2.5, -2.76, -2.64 and -2.60 were simulated for gamma-rays, protons, Neon 
and Fe respectively. All primary particles were assumed to be incident isotropically 
within the ~ 4°diameter field of view of the trigger pixels along a zenith direction 
of 50°. 

The large value of zenith angle was deliberately chosen to obtain a larger effective 
area of detection with a corresponding increase in the energy threshold as men- 
tioned earlier. Independent simulations of the TACTIC telescope carried out earlier 
have shown that the detection threshold for gamma-rays at > 45° zenith angle in- 
creases to ~ 10 TeV (2 TeV for vertical) while the effective collection area increases 
to ~ 10 9 cm 2 [12,17]. In this simulation study, we have also evaluated the effective 
collection area for 5TeV gamma-rays for the TACTIC telescope at zenith direction 
of 50° as ~ 7.07 x 10 s cm 2 . All the simulated Cerenkov images were subjected to 
similar cleaning and parameterization procedure in order to derive the Hillas im- 
age parameters [18] and the recently introduced fractal and wavelet parameters as 
described briefly below. 

3.1 Image Shape(Hillas parameterization 

The simulated Cerenkov images recorded at core distances varying from 10-400m 
have been parameterized following the well-known second order moment analy- 
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sis of the derived photo-electron counts in the image pixels. Hillas shape param- 
eters^], namely Length (L), Width(W), Azwidth (A), Distance (D), Miss (M), 
Orientation angle(o;) and size (S) were determined for each simulated image. The 
derived values were sorted as an independent set, for core distances < 400m for a 
comparative analysis. 



3.2 Fractal and Wavelet Parameters 

The longitudinal and lateral development of air showers, represented by the height 
of shower maximum, the number of relativistic electrons and muons, the angular 
and spatial distribution of Cerenkov photons and their total number varies with the 
primary mass. This produces differences in the structure of Cerenkov images which 
can be recognized through pattern recognition techniques. Fractal and wavelet pa- 
rameters quantify the Cerenkov image structure on different scale lengths, the for- 
mer characterizing the global structure of photon distribution in the image while 
the latter characterize more localized structures[19,20,24,25]. The basic motivation 
and the technique employed for fractal and wavelet analysis of Cerenkov images 
has been described in detail in[21-25]. We have followed an essentially similar pro- 
cedure to calculate fractal moments and dimensions as well as wavelet moments 
and dimensions of different order. 



3. 3 Fractal and Wavelet Analysis of Cerenkov Images 

There has been increasing interest in wavelet based approach in the recent years 
due to its success in several applications. The wavelet transform is a two parameter 
transform[26-28] . The two parameter domains of the wavelet tranform are dilation 
(scaling or compression) and translation (space). The advantage of using wavelets 
in the signal or structure analysis is the ability to perform local analysis of the signal 
or structure i.e to zoom on any interval of time and space. As such wavelet analysis 
has the ability to reveal the subtle details of the data which other analysis techniques 
fail to detect. Wavelet Analysis is regarded as the sequence of versatile filtering 
approach which will allow to examine the very minute stuctures left after filtering 
of the images at different scale lengths i.e to zoom at any interval of frequency ,time 
or space of the signal and enables us to detect some hidden aspect of the data or 
image and brings out corrections to the image analysis. 

The pattern recognition technique for gamma-hadron separation employs fractal 
and wavelet analysis of Cerenkov images to exploit differences in the structure 
of the recorded images due to differences in the longitudinal and lateral develop- 
ment of showers initiated by the two primary species. Showers initiated by pri- 
mary gamma-rays and protons produce different numbers of relativistic electrons 
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and muons and also have different heights (and lateral extents) of shower maxi- 
mum. These differences lead to structures in the Cerenkov image on different scale 
lengths which can be parameterized through multifractal and wavelet dimensions 
and moments of different order[ 19-24] . The usual procedure for calculating these 
parameters involves dividing the Cerenkov image into M = 4, 16, 64, 256 equally 
sized, non-overlapping cells and finding the number of CDC counts in each cell 
(6.5 CDC = 1 photo-electron). The multifractal moments of order q(q = 2— 6) are 
calculated from the following relation: 



Fractal moment G q (M) = ^ 



M rR , 



N 



(1) 



where N is the total number of CDC counts in the image, Kj is the CDC count in 
the ) th cell and q= 2, 3, 4, 5, 6 is the order of fractal moment. The wavelet moments 
are obtained as follows: 



Wavelet moment W q (M) = ^ 

3=1 



1 r A',., lyr" 



N 



(2) 



The fractal scale length is given by v = log 2 M 

Both the fractal moment G q as well as the wavelet moment W q exhibit a propor- 
tionality with the scale length which can be expressed as, 



G q (M) oc M Tq and W q {M) oc (3) 

The fractal and wavelet dimensions of order q are then determined from the follow- 
ing relations: 



Fractal dimension D q = — — (4) 

q - 1 

Wavelet dimension B„ = ^ 9 (5) 

q - 1 



where r q and (3 q are the slopes of the line between natural logarithm with scale 
length v of G q (M) and W q (M) respectively. Since wavelet dimensions are more 
sensitive to local structures in the photon distribution in Cerenkov images, they are 
expected to have better efficiency for gamma-hadron separation as proton-initiated 
showers are known to have a non-uniform structure due to contributions from pions 
and muons [20-25]. 



5 



Table 1 

Parameter domains and retention factors of different primary species for Cerenkov events 
at < 400m; A-Accepted, R-Rejected 







% 


% 


% 


% 


SNo 


Domain 


7 


P 


Ne 


Fe 


1 


0.06< L< 0.46 
0.06< W< 0.26 


91. OA 


95.0R 


97. OR 


99.0R 


2 


0.30< L< 0.58 
0.20< W< 0.48 


96.0R 


87.0A 


96.0R 


95. OR 


3 


0.40< L< 0.89 
0.40< W< 0.62 


99.5R 


80.0R 


48.0A 


78.5R 


4 


0.54< L<0.87 
0.45< W<0.69 


99.5R 


82.5R 


66.0R 


37.0A 



3.4 Results of Image Parameterization and Fractal-Wavelet Analysis 

Fig.l shows a composite plot of the frequency distribution of image shape parame- 
ters derived for showers initiated by 5-20 TeV gamma-rays, protons, Ne and Fe pri- 
maries with shower cores located at < 400m from the TACTIC imaging telescope. 
We find that the distributions for the shape parameters length and width (L,W) show 
a clear distinction between gamma-ray, proton and Ne/Fe induced showers, espe- 
cially a clear distinction between gamma-ray and hadron initiated showers. While 
on the other hand, the distributions for the Alpha(a), Distance(D) and Miss(M) pa- 
rameters do not exhibit any substantial dependence on the primary mass indicating 
that these parameters can have a limited application in primary mass segregation. 

This is expected because gamma-ray initiated Cerenkov images are known to be 
compact in shape as compared to the rather diffuse images produced by hadron 
initiated showers due to the presence of a large number of muons in the latter. 
For core distances < 400m, only Length(L), Width(W) parameters show the above 
mentioned difference while the other distributions tend to show no segregation. In 
Hillas parameterization as depicted in Table. 1, L and W parameter cuts show seg- 
regation between gammas and hadrons (p, Ne and Fe). Applying these cuts with 
appropriate values of L and W at different levels on Cerenkov images recorded by 
the telescope yields an acceptance for gamma-rays ~ 91% , for p ~ 87% , for Ne 
~ 48% and for Fe ~ 37% . We have tried various data cuts in order to find the 
optimum parameter domain which preferentially selects the maximum number of 
events initiated by one primary species alone. The results summarized in Table. 1 
show that, on the basis of data cuts involving L, W parameters, it is possible to 
broadly segregate the simulated events in terms of the primary mass. While this 
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Fig. 1. Image parameter distributions for gamma and hadron-initiated showers 

may not be useful as far as determining the absolute cosmic-ray mass composition 
is concerned, this procedure can be used to at least know if the cosmic-rays undergo 
enrichment by heavy nuclei as the energy increases beyond ~ 10 TeV. 
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As mentioned earlier the simulated Cerenkov images produced by the four primary 
species 7, p, Ne, Fe of different energies have also been subjected to analysis to 
calculate fractal and wavelet moments G q and W q as well as the multifractal and 
wavelet dimensions D q and B q of order q(-6 <q< 6)[17-20]. The peak value of 
fractal dimension Dq for the simulated Cerenkov images of 7, p, Ne, Fe lies at 
(0.4±0.02), (0.6±0.02), (0.7±0.03) and (0.8±0.03) respectively. Similarly Bq val- 
ues for these images lie at (5±0.6) for 7, (2.0±0.8) for p, (3.0±0.6) for Ne and 
(4 ±0.6) for Fe. The frequency distributions of D q and B q have been plotted in or- 
der to find any differences which could be exploited to segregate the events in terms 
of the primary mass. In this paper, we have considered the entire range of cosmic 
primaries i.e diffuse gamma-rays, protons, Ne and Fe. The gamma-rays considered 
in our simulation are within a field of view of ~ 4°. Fig.2 and Fig. 3 show a compos- 
ite representative frequency distribution of fractal and wavelet parameters namely, 
D 2 , Dq, B 2 and Bq. A careful visual examination of Fig.2 reveals that while the 
distributions for fractal parameters D 2 and Dq exhibit a clear difference between 
events of 7-ray and Fe origin with events initiated by p and Ne primaries falling in 
between, they do not promise high efficiency for segregating the events on the basis 
of their mass as the distributions are rather wide with significant overlaps among 
the four distributions due to large background contribution. 

It may be mentioned here that the same general features are exhibited by the dis- 
tributions for fractal dimensions of different order q, indicating that the fractal pa- 
rameters have poor efficiency for mass segregation of the Cerenkov events. On the 
contrary, the frequency distributions for the wavelet parameters are found to be 
very narrow and peaked. In one particular case they are distinctly separated from 
one another for the four primary species considered here. This is seen clearly in 
Fig. 3 where we have shown two representative distributions of wavelet dimensions 
B 2 and Bq. We find that the frequency distributions of Bq for 5 -20 TeV 7, p, Ne 
and Fe initiated events are completely segregated from one another with < 1% 
gamma- events having Bq =(5±0.6), < 1% Ne events having Bq =(3.0±0.6), > 
87% p events having Bq =(2.0±0.8) and < 1% Fe events with Bq =(4.0±0.6). The 
distributions for wavelet dimensions of other orders exhibit the same broad features 
as seen in the frequency distribution of B 2 indicating that, while these parameters 
can be used for 7- p segregation, they can not be employed for assigning a particu- 
lar primary mass to the events due to large overlapping of background events. As is 
clear from Table. 2, We find that the wavelet dimension Bq provides the maximum 
efficiency for segregating the Cerenkov events in terms of their mass on an event by 
event basis, making it possible to study cosmic-ray composition in the multi-TeV 
energy range with an imaging atmospheric Cerenkov telescope. 

The results so far obtained from simulations in each domain of Bq initiated by these 
four cosmic-ray primaries were used to estimate their flux values at different ener- 
gies. Various domains of wavelet dimensions B 2 and Bq for each species have been 
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Table 2 

Parameter domains and retention factors of different primary species for events at various 
Bq values; A-Accepted, R-Rejected 







% 


% 


% 


% 


S No 


Domain 


7 


P 


Ne 


Fe 


1 


1.6< B 6 < 2.4 


99.0R 


98.5A 


99.0R 


99.5R 


2 


2.8< B 6 < 3.6 


99.0R 


98.5R 


96.0A 


98.5R 


3 


3.8< B 6 < 4.6 


98.5R 


99.0R 


97.5R 


95. OA 


4 


5.2< B 6 <5.6 


91. OA 


98.0R 


97.0R 


96.0R 



identified from the simulated data. Out of 2,76,904 simulated events, 2,40,903 ±491 
events were retained in the proton domain with B 6 = (2.0±0.8); 2,492 ± 44 events 
are falling in Ne domain with B e =(3.0±0.6); 2,216 ± 41 events are recorded in 
Fe domain with B 6 = (4±0.6) and 155 ± 11 events are recorded in the 7 domain 
with B 6 = (5±0.6). These values lead to the flux estimation of (3.8±0.5)x 10~ 3 
m~ 2 s- 1 sr~\ (4.3±0.5)x 10~ 5 m~ 2 s" 1 sr~\ (3.9±0.5)xl0~ 5 m~ 2 s^sr" 1 and 
(2.9±0.4)xl0~ 6 m~ 2 s _1 sr _1 for protons, Ne, Fe and gamma-rays above lOTeV, 
15TeV, 20TeV, and 5TeV energy respectively.. 



4 Experimental Results 

The Cerenkov imaging telescope TACTIC has been used to search for TeV gamma- 
ray signals from several galactic and extragalactic sources and has successfully de- 
tected emission from the Crab Nebula, Mkn 421 and Mkn 501 [29,30,31]. The tele- 
scope is also used to monitor an off-source or background region where it tracks 
a region of sky away from the source direction in order to provide an independent 
measure of the cosmic-ray background rate and also to validate the gamma-ray sig- 
nal retrieval methodology. The fractal dimensions obtained from these Cerenkov 
images showed no potential in terms of segregation of events on the basis of mass 
because of significant overlaps of distributions due to large background contribu- 
tion. These simulated event values also yield a cosmic-ray proton composition of 
87% which is close to the experimentally measured proton abundance of ~ 89% at 
TeV energies. The gamma-ray to proton ratio determined from these values is also 
in close agreement with the value of 0.001 at TeV energies reported from other at- 
mospheric Cerenkov telescope systems[32] while as the protons out number other 
species by few orders of magnitude. 

We have also subjected Cerenkov images of various cosmic-ray primaries recorded 
by TACTIC during off-source observation runs to the above discussed fractal and 
wavelet analysis to investigate whether the B 6 parameter distribution leads to the 
expected cosmic-ray mass composition. Fig. 4 shows the frequency distribution of 
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Fig. 2. Fractal dimension distributions of Cerenkov images generated by gammas-rays and 
hadrons (p,Ne,Fe) 

parameter B 6 for both simulated and experimentally observed Cerenkov images 
recorded by the TACTIC telescope during a 37h long off-source exposure along 
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a direction making a zenith angle of ~ 50°. The total number of events recorded 
during Crab-off source scans was 2,76,900 corresponding to the event rate of ~ 
2.0 Hz, which translates into detection threshold energy of 4 TeV for cosmic -ray 
proton initiated events(~ 2 TeV for gamma-ray initiated events) as mentioned in 
the simulation process. While 2,35,635±485 events fall within the proton domain 
£ 6 (2.0±0.8); 1,788 ± 39 events are falling in Ne domain with B 6 = (3.0±0.6); 
1,538 ± 34 events are recorded in Fe domain with B e = (4±0.6) while as 127 
± 9 events fall in the gamma domain with £? 6 =(5±0.6). Taking into account the 
exposure time (37h), TACTIC effective area 7.07 x 10 8 cm 2 and the TACTIC field 
of view 3.8 x 10 _3 sr, the flux for various cosmic-primaries such as p, Ne, Fe and 
gammas above 5 TeV within 40% systematic error are found to be (3.7±0.5)x 
lO^m-V 1 sr~\ (3.0±0.5)x lO^m" 2 s" 1 sr~\ (2.7±0.4)x lO^m" 2 s^sr" 1 
and (2.4±0.4)x 10~ 6 m~ 2 s _1 sr _1 respectively [29]. The results obtained show good 
agreement between the simulated and experimental data. 



5 Summary and Conclusions 

A number of attempts have been made earlier to assess the potential of imaging 
atmospheric Cerenkov telescopes for separating different components of primary 
cosmic-rays in the energy region above ~ 10 TeV[32,33,34]. It has been shown 
that, while Hillas image parameters of atmospheric Cerenkov events recorded in an 
imaging Cerenkov telescope are able to separate protons from all other nuclei, the 
image parameters are rather inefficient in segregating nuclei of different masses. 
Our analysis confirms this fact even though we find that the image parameters like 
Width and Azwidth have the potential to separate the primaries into groups of light 
(protons and Ne) and heavy (Fe) nuclei. Several workers [34,35] have also investi- 
gated the efficiency of the temporal characteristics of atmospheric Cerenkov pulses 
for primary mass segregation. It has been reported that while some simple pulse 
shape parameters, like pulse rise time, full-width at half maximum and pulse de- 
cay time, carry information useful for cosmic-ray mass composition studies, the 
efficiency for primary mass segregation can not be quantified in view of the lim- 
ited simulation database. A detailed feasibility study of cosmic-ray composition 
measurements with Cerenkov telescopes using fractal parameters of the recorded 
Cerenkov images has been presented earlier[23]. 

It has been shown that the mass sensitivity of all the parameters increases with the 
image size, possibly because of the decreasing event by event fluctuations with in- 
creasing primary energy(~ size). However, a separation in more than two or three 
mass groups seems prohibited by the relatively poor mass resolution. Our results 
show that for primary hadrons of different masses with individual energies pre- 
selected to produce different amounts of Cerenkov light in the atmosphere, very 
efficient mass segregation is produced by the wavelet dimension B 6 . As already 
noted, we find that the frequency distributions of B 6 values for the three primary 



11 




Fig. 3. Wavelet dimension distributions of Cerenkov images generated by gammas-rays and 
hadrons (p,Ne,Fe) 

nuclear species (p, Ne and Fe), considered by us in the simulations, peaks at dis- 
tinctly separate values with negligible mutual overlap. To conclude, we have shown 
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Fig. 4. Wavelet dimension distribution of simulated and experimental data 

from our preliminary simulation-based study, that Cerenkov image shape parame- 
ters are able to separate primary cosmic -ray events into two or three groups of 
nuclei in terms of their mass but are inefficient in segregating the events in terms 
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of primary mass on an event by event basis. However, we find that the wavelet di- 
mension parameter B 6 of the Cerenkov image, has immense potential for event by 
event segregation(2 1,22,23) of the recorded events in terms of the primary mass. 
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